1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300000 iteration;
- 200000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.04 1.14
## LDevsum 1.02 1.08
## dh0 2.44 5.11
## dh1 1.74 3.37
## dh2 1.97 3.85
## dh3 1.42 2.48
## dl0 1.91 3.70
## dl1 1.13 1.50
## dl2 1.65 3.06
## dl3 1.44 2.33
##
## Multivariate psrf
##
## 1.87
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1231.45087 | 22092.3768 |
| DIC3 | 1159.66531 | 22131.0244 |
| PWAIC | 43.49456 | 240.7268 |
| WAIC | 1186.57378 | 22154.1588 |
2. Separate Modeling of High-Level
2.1 Trace plots for convergence check
The current MCMC setting is:
- 2000000 iteration;
- 100000 burn-in;
- 190 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## dh0 1.03 1.09
## dh1 1.04 1.04
## dh2 1.00 1.00
## dh3 1.03 1.05
##
## Multivariate psrf
##
## 1.02
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| H4 | |
|---|---|
| DIC | 1391.26658 |
| DIC3 | 1250.91902 |
| PWAIC | 92.64765 |
| WAIC | 1321.56593 |
3. Separate Modeling for Low-level
3.1 Trace plots for convergence check
The current MCMC setting is:
- 300000 iteration;
- 200000 burn-in;
- 100 thinning.
3.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## Devsum 1.00 1.01
## LDevsum 1.00 1.01
## dl0 1.02 1.07
## dl1 1.08 1.33
## dl2 1.01 1.03
## dl3 1.04 1.16
##
## Multivariate psrf
##
## 1.08
3.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
3.4 WAIC results
| L4 | |
|---|---|
| DIC | 22187.433 |
| DIC3 | 22235.386 |
| PWAIC | 295.985 |
| WAIC | 22267.482 |